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Abstract 


This  note  discusses  the  extension  of  the  Baum-Liu-Tesche  (BLT)  equation  to 
include  the  effects  of  propagating  EM  fields  and  their  coupling  to  transmission  line 
structures.  The  usual  BLT  equation  describes  the  propagation  of  voltage  (and  current) 
waves  along  transmission  lines.  In  the  extended  BLT  formalism,  both  voltage  and  E-field 
quantities  are  treated  together  and  this  permits  the  addition  of  scatterers  and  apertures 
in  the  models. 
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On  the  Addition  of  EM  Fieid  Propagation  and 
Coupiing  Effects  in  the  BLT  Equation 


1.  Introduction 

Recent  developments  in  electromagnetic  (EM)  field  sources  producing  high  power 
microwave  (HPM)  or  ultra  wide-hand  (UWB)  transient  signals  have  led  to  a  concern  about  the 
effects  of  these  environments  on  modern  digital  systems.  Such  excitations,  which  collectively 
are  referred  to  as  high  power  electromagnetic  (HPEM)  fields,  could  be  inadvertent,  like  the 
environment  produced  by  a  search  radar  illuminating  a  nearby  aircraft,  or  intentional,  such  as  a 
deliberate  attack  on  a  system  by  an  adversary  using  an  electromagnetic  weapon  [1].  In  either 
case,  the  effects  of  these  HPEM  fields  may  include  system  upset,  and  in  some  cases,  permanent 
damage. 

Eigure  1  illustrates  a  simple  example  of  an  electrical  system  excited  by  an  external 
HPEM  source.  This  energy  source  can  provide  either  a  narrow-band  pulsed  EM  field,  or  a  fast 
transient  field  that  couples  energy  onto  the  long  lines  attached  to  the  system.  Additionally,  the 
incident  field  excites  various  apertures  of  the  system  enclosure.  This  EM  energy  is  able  to 
penetrate  into  the  system,  is  distributed  within  the  internal  circuitry,  and  ultimately  arrives  at 
the  digital  components,  at  which  point  an  upset  or  failure  may  occur.  Of  course,  if  the  system  is 
well  shielded,  the  internal  components  may  not  be  affected  by  this  excitation;  however  all 
systems  must  have  imperfections  in  the  external  shielding  (for  access  hatches,  power 
penetrations,  etc.)  and  the  survival  of  such  a  system  is  never  100%  assured. 
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Figure  1.  Illustration  of  an  electrical  system  excited  by  an  external  EM  source,  causing 
upset  or  damage  in  the  internal  circuitry. 
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To  better  understand  the  possible  effeets  of  the  HPEM  environments  on  sueh  eleetrieal 
systems,  it  is  desired  to  develop  a  eomputational  modeling  approaeh  that  will  allow  an  analyst 
to  estimate  the  transient  responses  at  the  eomponent  level  -  both  in  amplitude  and  wave  form  - 
and  to  eonduet  trade  studies  to  examine  the  effeetiveness  of  various  system  design  ehanges  or 
shielding  eonfigurations.  The  purpose  of  this  report,  therefore,  is  to  outline  a  modeling 
approaeh  that  ean  be  used  for  this  purpose. 

The  HPEM  problem  posed  in  Eigure  1  is  similar  to  past  problems  involving  the 
proteetion  of  eleetronie  equipment  against  the  effeets  of  the  nuelear  eleetromagnetie  pulse 
(NEMP)  [2].  Although  the  NEMP  threat  is  lower  in  frequeney  (with  a  speetral  eontent  from 
several  kHz  to  about  100  MHz),  there  is  a  eommon  point  of  departure  for  the  analysis  that  ean 
be  used  in  the  present  ease.  In  the  NEMP  analysis,  the  modeling  of  eleetromagnetie  (EM)  field 
interaetion  with  a  eomplex  system  or  faeility  has  been  aeeomplished  by  the  introduetion  of  the 
eoneept  of  eleetromagnetie  topology.  This  idea  was  initially  developed  by  Baum  [3],  and  later 
formalized  in  the  literature  [4,  5,  6,  7].  This  methodology  permits  the  system  to  be  viewed  as  a 
eolleetion  of  eondueting  surfaees  that  attenuate  the  ineident  EM  field  environment  and 
effeetively  shield  the  interior  regions  from  the  external  field. 

Eor  the  sample  system  of  Eigure  1,  the  topologieal  representation  of  the  shielding 
surfaees  is  shown  in  Eigure  2.  These  shielding  surfaees  are  not  perfeet,  however,  as  there  may 
be  eonduetors  passing  through  the  surfaees,  and  external  EM  energy  ean  penetrate  into  the 
system  via  eurrent  and  eharge  flowing  on  the  eonduetors.  Moreover,  there  may  be  apertures 
(holes)  in  the  shielding  and  EM  fields  may  penetrate  into  the  interior.  Einally,  EM  field 
diffusion  through  the  shield  material  is  a  possibility,  if  the  shield  material  is  not  highly 
eondueting. 
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Figure!.  The  system  shielding  topology  diagram,  showing  the  EM  field  interaction 
paths  for  the  example  facility  of  Figure  1. 
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To  develop  a  computational  model  for  the  EM  responses  within  the  system,  the  concept 
of  the  interaction  sequence  diagram  can  be  introduced.  This  is  a  diagram  containing  the  main 
pathways  that  the  EM  energy  takes  as  it  penetrates  into  the  system.  Realizing  that  the  most 
important  paths  are  usually  the  “hard- wired”,  or  conductive,  paths \  a  considerable  amount  of 
work  has  been  conducted  on  developing  a  way  of  computing  the  behavior  of  current  and 
charge  on  interconnected  wire  networks.  As  an  example,  Eigure  3  illustrates  the  conducting 
penetration  paths  within  the  interaction  diagram  of  Eigure  2,  and  it  are  these  paths  that  the  flow 
of  EM  energy  must  be  modeled. 


External  EM  Environent 


Figures.  Simplified  EM  interaction  paths  for  the  sample  system,  as  obtained  by 
considering  only  the  conducting  penetration  paths. 

By  using  transmission  line  concepts,  which  result  from  a  simplification  of  a  complete 
solution  of  Maxwell’s  equations  [8],  a  computationally  simple  model  that  predicts  the 
propagation  of  EM  energy  into  the  system  has  been  developed.  This  model  is  based  on  the 
Baum-Eiu-Tesche  (BET)  equation,  which  is  a  matrix  equation  describing  the  behavior  of  the 
voltage  and/or  current  at  all  of  the  junctions  (or  interconnections)  of  the  conductors  in  the 
network  [6,  9,  10,  11].  In  a  sense,  this  equation  is  similar  to  a  node  analysis  performed  in 
conventional  circuit  theory,  except  that  in  the  BET  equation,  the  effects  of  EM  propagation 
along  the  individual  conductors  are  considered. 

Eigure  4  shows  a  simple  example  of  a  network  of  interconnected  transmission  lines  that 
can  be  analyzed  by  the  BET  equation.  In  this  case,  each  conductor  is  located  over  the  reference 


*  For  the  case  of  NEMP  modeling,  the  effects  of  the  EM  field  coupling  paths  in  the  interaction  sequence  diagram 
are  often  neglected.  This  results  in  considerations  of  EM  energy  flow  only  along  the  conductors  (e.g.,  the  hard¬ 
wired  lines).  Eor  the  higher  frequency  HPEM  analysis,  the  field  coupling  paths  may  be  more  important,  and  hence, 
they  must  be  included  in  a  more  general  model. 
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ground  plane  and  at  each  open  end  of  the  line  a  terminating  impedance  is  found.  This  network 
appears  as  if  it  were  constructed  by  single-wire  lines,  but  in  reality,  each  transmission  line  can 
be  a  multiconductor  bundle  (a  “tube”)  of  conductors  with  the  terminating  impedances  being 
generalized  N-port  loads.  This  diagram  illustrates  the  network  as  having  a  simple  tree-like 
configuration,  although  with  the  BLT  analysis  formalism,  it  is  also  possible  to  consider 
networks  having  one  or  more  loops  being  formed  by  the  transmission  line  tubes.  Details  of  this 
type  of  analysis  are  provided  in  ref.  [9]. 


Reference  Conductor 


Figure  4.  Example  of  a  transmission  line  network,  which  is  based  on  the  simplification 
of  the  interaction  sequence  diagram. 

Based  on  the  topological  modeling  and  the  BLT  equation,  a  computer  code  was  written 
for  analyzing  transmission  line  networks  like  that  of  Figure  4.  This  early  code  was  QV7TA  and 
its  operation  is  described  in  ref.  [12].  Subsequent  code  development  in  France  resulted  in 
another  code  -  CRIPTE  -  which  is  presently  a  commercially  available  code  for  performing  a 
wide -band  transmission  line  analysis.  The  details  of  this  code  and  typical  results  have  been 
described  by  Parmantier  [13]. 

Figure  5  illustrates  an  example  of  the  types  of  problems  that  can  be  analyzed  using  this 
modeling  technique.  Part  a  of  the  figure  shows  a  cable  harness  located  within  the 
Electromagnetic  Pulse  Test  Bed  Aircraft  (EMPTAC),  which  is  a  Boeing  707  aircraft  that  has 
been  converted  into  a  laboratory  to  validate  EMP  hardening  technologies.  This  aircraft  is 
located  at  the  Kirtland  APB  in  Albuquerque,  NM. 

As  a  result  of  a  collaboration  between  the  Office  National  d’Etudes  et  de  Recherches 
Aerospatiales  (ONERA)  (which  is  the  is  the  Prench  equivalent  of  NASA),  the  Centre  d’Etudes 
de  Gramat  (CEG)  and  the  Phillips  Eaboratory  of  the  US  Air  Porce,  a  combined  analysis  and 
measurement  program  was  conducted  to  examine  the  accuracy  of  the  CRIPTE  code.  Pigure  5b 
illustrates  the  network  model  developed  for  this  cable  harness,  and  a  comparison  of  the 
measured  and  computed  results  for  voltages  at  one  load  in  the  network  are  shown  in  part  c.  As 
can  be  noted,  the  agreement  between  the  model  and  measurements  is  quite  good. 
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shielded  cable 


b.  Transmission  line  network  for  the  aircraft  cable  harness. 


c.  Comparison  of  measured  and  calculated  cable  voltages,  shown  as  a  function  of  frequency. 

Figures.  Illustration  of  the  use  of  topological  modeling  concepts  for  cable  responses 
within  an  aircraft  (after  Parmantier  [13]). 
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Parmantier’s  work  on  the  modeling  of  the  aircraft  cable  harness  was  directed  towards  a 
NEMP  analysis  of  the  aircraft  equipment  response,  and  consequently,  neglecting  the  EM  field 
coupling  paths  is  a  reasonable  assumption.  However,  in  attempting  to  apply  these  topological 
models  to  higher  frequencies  that  are  typical  of  HPM  or  the  UWB  threats,  it  becomes  desirable 
to  include  these  additional  field-coupling  paths  in  the  BET  equation. 

To  show  what  is  desired  in  this  regard,  consider  the  diagram  in  Eigure  6,  which  shows  a 
simple  enclosure  with  an  external  antenna  with  a  conducting  penetration,  and  an  aperture 
penetration.  An  external  EM  source  produces  an  incident  field  that  couples  to  the  antenna, 
producing  currents  and  charges  that  penetrate  into  the  system  interior  by  propagation  along  the 
conductors  (the  solid  lines).  In  addition,  the  incident  field  is  able  to  penetrate  through  the 
aperture,  and  then  propagate  to  the  interior,  where  it  couples  to  an  internal  conductor.  This 
field  propagation  is  indicated  by  the  dotted  lines  in  the  figure.  Both  of  these  coupling  paths 
provide  a  response  at  an  internal  load  point. 

The  conventional  BET  equation  describes  the  propagation  of  EM  energy  (in  the  form  of 
voltage  or  current)  along  the  conductors  only,  with  the  excitation  arising  from  knowledge  of 
the  excitation  EM  field  acting  on  the  antenna.  A  more  general  form  of  the  BET  equation  can  be 
envisioned  that  models  the  propagation  of  EM  energy  both  on  the  conductors  and  through 
space,  with  the  excitation  function  now  being  a  specification  of  the  fundamental  radiation 
properties  of  the  external  EM  source  -  not  just  an  incident  EM  field  on  the  system.  In  this 
way,  a  more  complete  EM  interaction  model  results,  and  its  responses  are  expected  to  be  more 
accurate  at  the  higher  frequencies. 


Antenna 


Figure  6.  A  simple  example  of  a  more  general  interaction  path,  containing  both  EM  field 
and  hard-wired  interaction  paths. 

This  report  serves  to  outline  the  approach  to  be  taken  to  extend  the  BET  equation  to 
account  for  the  effects  of  EM  field  propagation  within  the  system.  We  take  a  top-level 
approach  to  show  the  general  extension  to  the  model,  and  in  doing  so,  we  neglect  a  number  of 
effects  that  can  be  of  importance  in  the  application  of  this  model  to  real  systems.  Specifically, 
we 
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•  Neglect  the  near-field  effects  in  the  behavior  of  the  source  and  the  transmission  lines. 
In  doing  so,  only  the  far-field  plane-wave  relationships  for  the  coupling  and  scattering 
of  the  elements  in  the  model  are  considered. 

•  Neglect  the  possibility  of  different  polarizations  of  the  EM  fields 

•  Neglect  the  possibility  of  the  existence  of  a  common  mode  (or  “antenna”  mode) 
response  on  the  transmission  lines.  This  implies  that  only  the  transmission  line  mode  is 
considered,  which  is  suitable  for  the  case  when  the  conductor  is  located  over  a 
reference  groundplane. 

•  Neglect  the  higher  order  modes  that  can  exist  on  the  conductors,  and  assume  that  the 
coupled  EM  field  responses  (the  line  current  and  charge)  are  adequately  represented  by 
transmission  line  theory. 

•  Assume  that  the  scattered  EM  field  from  any  of  the  conductors  within  the  network  is 
adequately  modeled  from  knowledge  of  the  currents  obtained  from  transmission  line 
theory. 

What  we  concentrate  on  in  this  paper  is  the  manner  in  which  the  BET  equation  must  be 
modified  to  take  into  account  the  field  coupling.  While  including  all  of  the  above  effects  would 
certainly  be  possible,  we  believe  that  they  might  obscure  the  main  features  of  the  BET 
extensions,  and  consequently,  these  refinements  are  postponed  for  future  discussion. 

The  starting  point  of  the  present  paper  is  a  short  review  of  the  derivation  of  the  BET 
equation  for  a  simple  section  transmission  line.  This  is  done  in  Section  2,  which  provides  an 
introduction  to  the  BET  concept  and  shows  how  it  can  be  applied  to  a  more  extensive  network 
like  that  of  Eigure  4.  Then,  in  Section  3  the  steps  necessary  to  extend  the  BET  equation  to 
include  the  field  coupling  paths  is  described.  This  methodology  is  illustrated  with  a  very 
simple  example  of  a  single  dipole  source  coupling  to  a  transmission  line.  In  Section  4,  the  BET 
theory  is  further  extended  to  treat  the  problem  of  field  coupling  from  a  dipole  source  through 
an  aperture  to  a  transmission  line  and  a  simple  example  is  given.  Einally,  in  Section  5  these 
results  are  summarized  and  references  are  provided  in  Section  6.  Appendices  A  and  B  contain 
further  information  on  the  field  coupling  and  radiation  functions  for  the  transmission  line 
structure  used  in  the  examples. 

It  should  be  noted  that  there  are  no  numerical  calculations  or  illustrations  of  this  theory 
in  the  present  paper.  The  goal  of  the  discussion  here  is  to  outline  the  analysis  methodology, 
and  in  doing  so,  many  of  the  details  of  field  polarization  and  spatial  variations  of  the  radiated 
or  scattered  fields  are  neglected.  We  anticipate  that  these  details  will  be  filled  in  at  a  later  date 
in  a  companion  paper  that  provides  numerical  demonstrations  of  this  BET  formalism. 
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2.  Review  of  the  Derivation  of  the  BLT  Equation 

To  begin  the  development  of  the  BLT  equation  containing  distributed  EM  field  signal 
paths  and  excitations,  let  us  first  review  the  basic  steps  used  in  deriving  the  conventional  BLT 
equation  for  transmission  line  networks.  Instead  of  treating  a  general  transmission  line 
network,  however,  we  outline  the  formulation  of  the  equation  for  a  single  transmission  line,  as 
shown  in  Figure  7.  This  development  will  shown  the  overall  analysis  approach,  which  can  be 
generalized  for  a  network  of  interconnected  multiconductor  transmission  lines,  as  detailed  in 
refs.  [2]  and  [9].  This  analysis  is  conducted  in  the  frequency  domain,  with  an  assumed  time 
dependence  of  e  Recall  also  that  we  are  limiting  our  discussion  to  transmission  line 
currents  only:  the  antenna  mode  is  not  considered  here.  Physically  this  implies  that  we  are 
considering  a  single  conductor  line  located  over  a  groundplane,  with  the  image  conductor  in 
the  groundplane  serving  as  the  second  conductor  in  the  2-wire  line. 

The  geometry  shown  in  Figure  7  is  that  of  a  simple  uniform  two- wire  transmission  line 
of  length  L  and  characteristic  impedance  Z^.  We  assume  that  it  has  a  propagation  constant  % 
which  in  the  most  general  case,  can  be  dispersive  if  the  line  has  resistive  losses.  The  line  is 
terminated  at  each  end  by  impedance  loads  ^L1  and  Z/^2’  it  is  assumed  to  be  excited  by  a 

lumped  voltage  and/or  current  source,  Vg  and  1^,  at  the  position  v  =  along  the  line. 


General  expressions  for  the  induced  current  I(x)  and  voltage  V(x)  at  any  position  x 
along  the  line  have  been  developed  in  [8].  In  many  instances,  however,  only  the  load  voltages 
Vj  and  V2,  or  the  corresponding  load  currents,  are  desired  for  the  analysis.  In  this  case,  the 
results  of  performing  an  analysis  for  these  load  quantities  can  be  put  into  a  very  compact  form, 
known  as  the  BLT  equation. 


Transmission 

Line 


Node  #1 


Excitation 


Observation 

point 


Node  #2 


Linear 

Graph 


Node  #1 


Forward  traveiing 
wave 

vqx) 


(X) 

Backward  traveiing 
wave 


Excitation 


Node  #2 


Figure  7.  A  simple  transmission  line  with  lumped  source  excitations  (top)  and  its  linear 
graph  representation  (bottom). 
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2.1  The  Voltage  BLT  Equation 

While  a  BLT  equation  can  be  developed  for  either  the  load  voltages  or  the  load 
currents,  this  development  will  concentrate  on  the  equation  for  the  voltage  responses.  The 
starting  point  of  this  analysis  is  to  recognize  that  there  are  voltage  traveling  waves  propagating 
in  the  positive  direction  (y+)  and  in  the  negative  direction  (V"),  as  shown  in  the  bottom  portion 
of  Figure  7.  Application  of  transmission  line  theory  to  this  problem  shows  that  these  waves 
propagate  on  the  line  with  a  functional  dependence  of  e^^ ,  respectively,  and  at  any  location  x 
along  the  line,  the  total  line  voltage  is  calculated  as  the  sum  of  the  traveling  waves,  V{x)  = 
y+(v)  +  V^{x).  The  term  yis  the  propagation  constant  of  the  waves  along  the  transmission  line, 
which  to  first  order  is  equal  to  jk  =  jlnflc  ,  which  is  the  propagation  constant  of  waves  in  free- 
space. 


As  noted  in  Figure  7,  it  is  also  useful  to  consider  describing  the  voltage  waves  at  the 
junctions  by  the  incident  and  reflected  voltage  waves,  and  Note  that  at  nodes  1  and  2 
of  the  simple  transmission  line  graph,  we  have  the  following  definitions: 


y^(0)  =  yi^ 

;  V-{Q)^Vf=Vr 

(at  node  1) 

(la) 

V^{L)=Vf  =vr  : 

;  V-{L)^Vf=Vff 

(at  node  2) 

(lb) 

For  this  simple  transmission  line,  the  voltage  waves  are  initially  excited  by  the  lumped 
voltage  and  current  source  at  x  =  x^.  As  suggested  in  Figure  8  the  sources  produce  a  positive 
traveling  wave  for  x>  x^  and  a  negative  traveling  wave  for  v  <  Ref.  [8]  derives  expressions 
for  these  voltage  waves  in  terms  of  the  source  strengths  as 

=  =  0  (forx>vJ  (2a) 

=  ;y;(.^)  =  0  (forx<vJ  (2b) 


Node  #1  Excitation  Node  #2 


Figure  8.  Illustration  of  the  positive  and  negative  traveling  waves  excited  by  the  lumped 
sources  at  jc  = 
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With  no  excitation  sources,  the  positive  traveling  voltage  at  junction  #2  is  expressed  in 
terms  of  the  positive  traveling  wave  at  junction  #1  as 

(3) 

However,  with  the  source  present,  there  is  an  additional  component  to  that  is  due  to  the 
positive  wave  traveling  away  from  the  source.  In  this  manner,  Eq.(3)  is  modified  to  include 
the  source  terms,  as 


y;  =y,^e-^  +i(y^  (4) 

Similarly,  the  negative  traveling  wave  at  junction  #1  may  be  expressed  in  terms  of  the 
negative  wave  at  junction  #2,  plus  the  source  contribution  as 

(5) 

Expressing  the  voltage  traveling  waves  in  terms  of  the  incident  and  reflected  waves  at 
the  junctions  using  Eq.(l),  Eqs.(4)  and  (5)  can  be  put  into  a  matrix  equation  for  the  incident 
waves  at  the  two  junctions  as 


inc 

1 - 

1 

o 

1 _ _ 

_ 1 

yinc 

- 1 

o 

1 

_ _ _ 1 

1 

■ 

_ 1 

+ 

Rearranging  the  terms  of  Eq.(6)  provides  a  similar  expression  for  the  reflected  voltages  at  the 
loads: 


yref 
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O 

_ 1 

inc 

1 

1 
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+ 

N 

_ 1 

- 1 

o 

_ 1 

T  r  inc 
-2 

1 

1 

N 

1 

-|r.| 

+ 

_ 1 

(7) 


Transmission  line  theory  [8]  provides  a  useful  relationship  between  the  incident  and 
reflected  voltage  waves  at  a  load  through  the  expression 

(g) 

where  p  is  the  voltage  reflection  coefficient.  Eor  an  impedance  load  on  a  transmission  line 
having  a  characteristic  impedance  Z^,  the  reflection  coefficient  is  defined  as 


^  In  mathematical  terms,  this  expression  is  viewed  as  a  solution  to  the  wave  equation  with  sources,  which  contains 
a  homogeneous  solution  and  a  particular  solution. 
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(9) 


For  this  simple  transmission  line  network,  the  vector  of  reflected  voltages  in  Eq.(7)  can 
be  expressed  in  terms  of  the  vector  of  incident  voltages  by  the  matrix  relationship 


1 

1 

>1 

0  ■ 

p7  inc 

1 

1 

0 

Pi. 

(10) 


where  pi  and  p2  are  the  reflection  coefficients  at  each  end  (node)  of  the  line. 


Substituting  Eq.(lO)  into  Eq.(7)  and  solving  for  the  vector  of  incident  voltages  provides 
the  first  form  of  the  voltage  BLT  equation: 
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_ 1 
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or  upon  combining  the  terms  of  the  inverted  matrix, 
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(11) 


(12) 


This  is  an  equation  for  the  incident  voltage  waves  at  both  ends  of  the  line,  in  terms  of  the 
transmission  line  parameters,  the  load  impedances  and  the  source  excitations. 


The  second  form  of  the  voltage  BLT  equation  is  for  the  total  load  voltage  vector,  which 
is  obtained  using  Eq.  (10),  as 


1^1 

.^2 
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(13) 


Substituting  this  expression  into  Eq.(12)  yields  the  following  total  voltage  BLT  equation: 


'^1  ' 

'1  +  A 

0 

-A 

-1 

■^r 

.^2: 

0 

1  + A_ 

-Pi. 
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where  the  source  (excitation)  term  for  the  lumped  current  and/or  voltage  sources  at  v  =  is 
defined  as 


1 

_  1 
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-1- 

N 

_ 1 
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1 _ 

1 
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N 
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-|(N 
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_ 1 

(15) 


Note  that  Sj  represents  the  source  term  for  the  excitation  of  the  positive  traveling  voltage  wave 
on  the  line,  while  S2  denotes  the  source  excitation  for  the  negative  traveling  voltage  wave. 


2.2  The  Current  BLT  Equation 

A  similar  analysis  involving  current  traveling  waves  can  be  used  to  derive  a  BLT 
equation  for  the  currents  flowing  in  the  load  impedances  of  the  line.  In  doing  so,  care  must  be 
used  in  accounting  for  the  difference  in  sign  between  the  line  current  and  the  load  current  at 
the  load  at  v  =  0.  Alternatively,  one  can  use  the  load  voltage  BLT  equation  (14)  and  the  V-I 
relationship  at  the  loads  to  determine  the  load  current.  In  any  case,  the  following  BLT 
equation  for  the  loads  currents  results: 


a' 

1 

'l-A 

0 

.  “A 

e^  " 

-1 

h_ 

0 

1-A_ 

-P2_ 

(16) 


2.3  The  BLT  Equations  for  Incident  Field  Excitation 

The  BLT  equations  for  the  load  voltages  and  currents  due  to  lumped  sources  at  x  = 
derived  in  the  previous  sections  can  be  used  to  determine  the  induced  load  responses  arising 
from  the  illumination  of  the  line  by  an  incident  EM  field,  as  shown  in  Figure  9.  As  noted  in  [7] 
and  [14],  an  incident  EM  field  acting  on  a  two-wire  line  can  be  thought  of  as  producing  voltage 
sources  that  are  distributed  along  both  of  the  wires  of  the  line,  as  well  as  two  lumped  voltage 
sources  at  the  line  ends  .  The  distributed  voltage  sources  induce  both  the  antenna  mode  (which 
we  neglect  in  this  discussion)  and  the  transmission  line  mode,  which  contributes  to  the  load 
currents. 


^  This  analysis  is  based  on  the  so-called  Agrawal  formulation  of  the  transmission  line  excitation  and  is  formally 
valid  for  a  transmission  line  in  free  space  or  in  a  homogeneous  dielectric  material.  For  the  interesting  case  of  a 
dielectric  coated  transmission  line,  the  use  of  the  Agrawal  formulation  is  open  to  question,  and  is  a  current 
problem  of  interest  in  transmission  line  theory. 
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Figure  9.  Field  excitation  of  a  two-wire  transmission  line,  showing  the  voltage  excitation 
sources  needed  in  the  BLT  equation. 

The  derivation  of  the  field  coupling  to  the  simple  two-wire  line  is  presented  in 
Appendix  A,  where  it  is  shown  that  the  same  BLT  equations  (14)  and  (16)  for  the  load  voltages 
and  currents  can  be  used  to  describe  the  field  coupling  responses,  but  with  a  different  source 
vector.  Assuming  that  the  propagation  constant  on  the  line  is  the  same  as  that  of  free  space  (/= 
jk  =  jcc/c),  this  source  vector  is  found  by  integrating  the  effects  of  the  distributed  field  sources 
along  the  line  and  combining  the  effects  of  the  end  sources. 

For  the  special  case  of  the  incident  field  being  a  plane  wave  with  an  angle  of  incidence 
\\f  and  the  incident  E-field  lying  in  the  plane  of  the  transmission  line,  the  resulting  source  term 
has  the  following  form"^: 


_  E‘’“^d 

i^ikm-cosy/)  _|j 

2 

^jkL^-jkL{\+cmy/)  _ 

=  E 


inc 


(17) 


Inserting  Eq.(17)  into  Eq.(14)  or  (16)  thus  provides  the  load  responses  for  plane  wave 
excitation  of  the  line. 


2.4  Generalized  Forms  for  the  BLT  Equations 

In  examining  the  derivation  of  the  BET  equation  for  the  single  transmission  line,  it  is 
possible  to  illustrate  the  forms  that  the  solutions  take  for  more  general  networks  or  for 


Ref. [8]  presents  the  more  general  expression  for  this  source  vector  for  the  case  of  an  arbitrarily  incident  and 
polarized  EM  field  acting  on  the  line,  and  for  the  case  of  the  transmission  line  propagation  constant  ^differing 
from  the  free  space  propagation  constant,  ']k.  This  results  in  a  more  complicated  expression  for  the  excitation 
terms. 
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multiconductor  transmission  lines.  Consider  the  definitions  of  the  following  terms,  which  have 
been  used  in  the  previous  sections: 


Line  Propagation  Matrix  F  = 


Load  Scattering  Matrix  S  = 


0 


Characteristic  Admittance  Matrix  F 


Unit  Matrix  U  = 


0 


A  0 
0  P2 


1/Z,  0 

0  1/Z„ 


1  0 
0  1 


7, 

Response  Vectors 

V  = 

;  /  = 

1 

^2^ 

_^2_ 

(18) 

(19) 

(20) 

(21) 

(22) 


With  these  definitions,  the  total  voltage  BLT  equation  (14)  takes  the  shorthand  form 


F 


U  +5 


'r-5f 


and  the  corresponding  current  BLT  equation  is 


/  =  F 


U -S 


r-s 


'-1  . 


(23) 


(24) 


As  discussed  in  [9]  the  BLT  equation  for  a  network  of  interconnected  transmission 

lines  has  the  same  form  as  Eq.(23)  or  (24).  The  matrix  F  represents  all  of  the  propagation 
functions  (in  the  forward  and  backward  directions)  on  each  of  the  transmission  lines  in  the 
network.  Thus,  if  there  are  N  transmission  lines  (sometimes  called  “tubes”)  in  the  network, 

then  F  is  of  order  2N  x  2N.  Assuming  that  the  responses  are  ordered  according  to  the  tube 

number  in  the  network,  F  is  a  block  diagonal  matrix,  with  each  sub-matrix  defining  the 
propagation  terms  for  the  individual  transmission  lines. 

The  scattering  matrix  S  contains  all  of  the  scattering  coefficients  from  all  of  the 
junctions  in  the  network.  For  the  responses  ordered  by  the  tube  number,  this  matrix  is  sparse, 
but  not  necessarily  block  diagonal,  since  the  locations  of  the  various  scattering  coefficients 
depend  on  how  the  junctions  in  the  network  are  numbered  and  interconnected. 
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3.  Extension  of  the  BLT  Equation  for  Distributed  Field  Excitation 

In  goal  of  this  paper  is  to  extend  the  eoneept  of  the  transmission  line  network  and  its 
solution  via  the  BLT  equation  to  the  distributed  field  problem  shown  in  Figure  10.  This 
configuration  consists  of  a  source  of  EM  energy  located  outside  a  shielded  enclosure,  which 
contains  an  internal  transmission  line.  Generally,  in  the  external  and  internal  regions  there  can 
be  other  scattering  objects  that  influence  the  fields.  The  shield  is  assumed  to  have  a  small 
aperture;  the  EM  fields  that  leak  through  this  aperture  are  able  to  induce  currents  and  charge  on 
the  wire  and  the  load  voltage  and  current  responses  are  to  be  determined.  Reference  [15]  has 
discussed  one  aspect  of  this  problem,  and  here  we  will  develop  an  alternate  approach  that  is 
consistent  with  the  assumptions  of  transmission  line  theory. 


Excitation  Region 
"outside" 


EM  Field 
Source(s) 

P&M 


Shielded  Region 

"inside"  i-'W-i  Load  #1 


Induced  Load 
Responses 


Scattering 

object 


Load  #2 


EM  Field 
Source 


Shield  and 
Aperture 


Transmission  Line 


O 

Scattering 

object 


Figure  10.  Problem  geometry  for  an  external  EM  field  source  producing  a  current 
response  on  a  transmission  line  witbin  a  shielded  region. 


As  mentioned  earlier,  we  are  not  interested  in  providing  a  detailed  numerical  solution  at 
this  stage.  Consequently,  we  will  ignore  many  of  the  (important)  subtleties  of  such  an  EM 
analysis,  including  EM  field  polarization  and  angle  of  incidence,  geometrical  details  of  the 
shield  and  aperture,  source  details  (radiation  pattern,  etc.),  and  the  fact  that  antenna  modes  can 
be  induced  on  the  transmission  line  structure  in  addition  to  the  transmission  line  modes.  We 
wish  to  avoid  these  details  and  concentrate  on  the  basic  concepts  of  extending  the  BET 
equation  to  be  able  to  treat  EM  field  propagation  and  coupling. 

Eollowing  the  development  of  the  BET  equation,  as  summarized  in  Section  2,  we  wish 
to  describe  the  EM  interaction  problem  of  Eigure  10  using  a  description  of  traveling  waves  on 
appropriate  tubes  (or  branches)  and  reflection  coefficients  at  the  appropriate  junction  (or 
nodes)  within  the  problem  volume.  As  is  evident  from  this  figure,  in  some  parts  of  the 
problem,  the  traveling  waves  can  describe  electromagnetic  fields  that  propagate  from  one  point 
in  space  to  another.  Traveling  waves  also  exist  on  the  transmission  line,  and  in  this  case,  they 
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are  best  deseribed  by  voltage  and  eurrent  disturbanees,  rather  than  by  EM  fields.  These  latter 
traveling  waves  are  those  already  treated  by  the  eonventional  BET  equation. 

Similarly,  at  some  loeations  in  the  problem,  refleetions  ean  be  represented  by  EM  field 
refleetion  eoeffieients,  and  at  other  loeations,  transmission  line  voltage  or  eurrent  refleetion 
eoeffieients  are  appropriate. 

A  fundamental  issue  in  developing  a  BET  representation  for  the  interaetion  problem  of 
Eigure  10  is  how  to  treat  the  EM  field  eoupling  between  the  ineident  field  and  the  transmission 
line.  To  this  end,  it  is  prudent  to  solve  a  simple,  preliminary  problem  that  shows  how  this 
eoupling  ean  be  ineorporated  into  the  BET  equation,  and  then  later  eonsider  the  solution  to  the 
more  eomplex  problem  of  Eigure  10. 

To  this  end,  the  simple  souree/transmission  line  problem  shown  in  Eigure  11  ean  be 
eonsidered.  A  simple  transmission  line  having  load  speeified  at  eaeh  end  is  illuminated  by  a 
distant  eleetrieal  souree,  whieh  induees  eurrents  and  voltages  at  the  loads.  In  this  simple 
example,  we  negleet  the  presenee  of  the  seattering  objeets  shown  in  Eigure  10,  and  eoneentrate 
only  on  the  exeitation  of  the  transmission  line  and  its  re -radiation. 

By  knowing  the  ineident  E- field  at  the  transmission  line,  the  BET  equations  (14) 

or  (16),  together  with  the  exeitation  terms  in  Eq.(17)  ean  be  used  direetly  to  determine  the  load 
responses  (for  assumed  plane  wave  exeitation).  In  this  manner,  the  ineident  E- field  is 
eonsidered  to  be  the  primary  exeitation  of  the  problem  (as  this  exeitation  oeeurs  in  the 
“foreing”  term  of  the  BET  equation. 


Figure  11.  Geometry  of  a  simple  elementary  source  exciting  a  transmission  line. 
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What  we  desire,  however,  is  to  consider  the  excitation  in  the  problem  shown  in  Figure 
11  to  be  the  dipole  sources  themselves  [P  and  m),  not  the  incident  fields  they  produce.  In  this 
way,  the  incident  field  on  the  transmission  line,  along  with  a  scattered  EM  field  from  the  line 
will  constitute  additional  response  variables  in  an  expanded  BET  equation. 

The  analysis  of  this  simple  problem  can  be  described  in  terms  of  the  signal  flow 
diagram  illustrated  in  Eigure  12,  which  shows  the  paths  of  energy  flow.  The  transmission  line 
is  represented  by  the  “conventional”  tube  with  solid  lines  and  is  denoted  as  Tube  1.  This  is 
identical  to  the  usual  tubes  used  in  the  BET  equation,  with  forward  and  backward  propagating 
voltage  and  current  waves.  Nodes  1  and  2  represent  the  transmission  line  loads,  as  described  in 
Section  2,  and  these  are  represented  by  reflection  coefficients  relating  the  incident  and 
reflected  voltage  and  current  waves  at  these  locations. 

In  addition  to  the  transmission  line  tube,  there  is  an  EM  field  propagation  “tube”,  which 
is  denoted  by  the  dotted  lines.  This  is  Tube  2  and  it  represents  the  path  by  which  the  EM  field 
propagates  from  the  source  to  the  line  (at  node  4),  and  also  to  a  distant  field  observation  point, 
at  node  3.  At  these  nodes,  the  responses  that  are  to  be  included  in  the  BET  equation  are  the  E 
and  H  fields. 

Nodes  1,  2  and  3  can  be  thought  of  as  conventional  nodes,  in  which  there  are  incident 
and  reflected  voltages  and  currents  (or  E  and  H-fields).  The  field  node  3  is  really  not  of  interest 
here,  because  it  might  be  thought  of  as  being  “at  infinity”  with  a  zero  reflection  coefficient  to 
represent  outward  only  propagation  of  the  EM  fields  in  free  space.  However,  in  more 
complicated  problems,  this  node  may  be  important,  say  if  there  were  an  electrical  conducting 
wall  or  other  object  far  from  the  source  and  transmission  line.  Thus,  this  node  will  remain  in 
the  sample  problem  discussed  here. 

Node  4  is  a  different  type  of  node  from  the  others,  in  that  it  contains  the  EM  field 
coupling  and  scattering  information  for  this  problem.  Consequently,  it  is  represented  by  a 
different  type  of  node  symbol  -  a  box  rather  than  a  circle.  At  this  node,  there  is  an  incident  EM 
field  component  produced  by  the  dipole  sources,  as  well  as  a  reflected  EM  field  from  node  3. 
This  incident  EM  field  is  “reflected”  by  the  transmission  line  structure,  and  results  in  a 
scattered  EM  field  propagating  away  towards  node  3.  At  the  same  time,  this  incident  EM  field 
on  the  line  contributes  to  the  incident  voltage  and  current  waves  at  the  transmission  line  nodes 
1  and  2.  Eurthermore,  the  reflected  voltage  and  current  waves  from  nodes  1  and  2  are  able  to 
re-radiate,  and  this  produces  a  component  of  the  scattered  EM  field  from  the  transmission  line. 
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Node  2 


Figure  12.  The  signal  flow  graph  for  the  EM  interaction  geometry  of  Figure  11. 


To  begin  with  the  development  of  the  modified  BLT  equation,  we  will  first  clarify  a 
notational  issue.  A  traveling  wave  component  (voltage,  current,  E- field,  or  H-field)  will  be 
found  on  each  tube  i  and  at  each  node  j  of  the  network.  This  will  be  represented  using  the 
notation  j.  Moreover,  for  this  problem,  we  will  concentrate  on  voltage  and  E-field  responses 
(with  the  corresponding  current  and  H-field  responses  being  obtained  from  characteristic 
impedances  of  the  transmission  line  and  free  space). 

To  develop  the  network  BET  equation,  we  will  need  to  examine  both  the  tube 
propagation  relationships  and  the  node  reflection  relationships,  as  described  in  the  next 
section. 


3.1  Tube  Propagation  Relationships 

Considering  first  the  transmission  line  Tube  1,  the  relationship  between  the  reflected 
and  incident  voltages  at  nodes  1  and  2  can  be  expressed  from  Eq.(7)  as 


1 

"  0 

T/  inc 
*^1,1 

-1- 

yref 

_ 

0 

_-f,{Ys)ez_ 

where  the  effects  of  the  incident  E-field  at  the  wire  ( )  can  be  represented  by  the  integral 
functions  F^iiPs)  F2(^J developed  in  Appendix  A.  Note  that  these  functions  depend  on 
the  local  angle  of  incidence  at  the  line  and  in  the  more  general  case,  on  the  polarization 
angles.  Note  that  Eq.(25)  contains  two  terms.  The  first  is  a  propagation  relationship  between 
incident  and  reflected  voltage  waves  on  the  transmission  line  having  a  propagation  constant  y. 
The  second  is  a  term  that  relates  to  the  excitation  of  waves  on  the  line  due  to  the  incident  field 
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A  similar  equation  for  incident  and  reflected  E-fields  in  tube  2  can  be  developed. 
Neglecting  polarization  issues,  the  incident  field  at  node  4  can  be  expressed  as 
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where  the  terms  take  into  account  the  spherical  propagation  of  fields  through  space  and 
denotes  the  amplitude  and  angular  source  term  from  the  elementary  dipole  sources, 
which  are  the  primary  sources  in  this  example.  The  term  y/j  represents  the  angle  of  departure 
of  the  radiation  from  the  source  to  the  transmission  line,  relative  to  some  convenient  source 
axis,  as  noted  in  Figure  11.  Furthermore,  the  term  aj  in  Eq.(26)  represents  a  typical  dimension 
of  the  scattering  body  at  node  3. 


Similarly,  at  node  3,  the  propagation  relationship  for  the  incident  E-field  can  be 
expressed  as 
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The  source  term  S2(^^2)  this  expression  is  the  same  as  in  Eq.(26),  but  with  a  different 
intensity  because  the  observation  angle  y/2  and  distance  from  the  source  to  node  3  is  different 
than  from  node  4.  The  incident  field  at  node  3,  »  ^l^o  contains  source  components  from 

the  traveling  waves  existing  on  the  transmission  line,  as  developed  in  Appendix  A,  by  way  of 
the  F j  and  F2  functions.  The  term  denotes  the  overall  size  of  the  node  4,  which  is  needed 
for  dimensional  consistency. 


Eqs.(26)  and  (27)  can  be  put  into  matrix  form  as 
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and  upon  rearranging  terms,  this  equation  can  be  solved  for  the  reflected  fields  at  the  nodes  as 
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Combining  the  propagation  portions  of  Eq.(25)  and  Eq.(29)  into  a  single  4x4  matrix 
for  the  reflected  quantities  provides  the  following  propagation  equation  for  this  composite 
problem: 
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(30a) 


Notice  that  the  equation  contains  the  information  about  how  voltage  wave  propagate  on 
the  transmission  line,  how  the  E-fields  propagate  in  free  space,  how  the  fields  couple  to  the 
transmission  line,  and  how  the  transmission  line  currents  produce  a  radiated  field. 


Eq.(30a)  is  seen  to  have  incident  and  reflected  variables  that  are  of  mixed  units:  volts 
and  volts/meter.  This  equation  can  be  put  into  a  more  consistent  form  involving  only  the  units 
of  volts,  by  multiplying  the  last  two  rows  of  the  equation  by  aj  and  a^,  respectively.  This 
yields  the  following  expression: 
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3.2  Node  Reflection  Relationships 

Eollowing  the  procedure  used  for  the  simple  transmission  line  in  Section  2,  we  now 
develop  an  expression  for  the  reflected  waves  using  the  node  parameters.  At  nodes  1  and  2,  the 
relationships  between  the  incident  and  reflected  voltages  are  given  by  the  reflection 
coefficients  on  the  transmission  line,  involving  the  characteristic  impedance  of  the  line  and 
the  load  impedance  Z^,  as 
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with 


23 


(32) 


where  now  relates  incident  and  reflected  E-fields  at  node  3.  As  previously  mentioned,  we 
assume  that  this  node  is  at  infinity  and  there  is  no  field  reflection  at  this  point,  so  =  0. 

The  definition  of  the  reflection  coefficient  at  node  4,  P4,  is  less  obvious  than  at  node  3. 
At  node  4,  it  is  clear  that  an  incident  field  at  the  node  will  be  reflected  by  the  transmission  line. 
However,  this  effect  is  already  accounted  for  in  the  propagation  relationships  of  Eq.(30).  The 
reflection  coefficient  at  node  4,  therefore,  is  the  reflection  coefficient  in  the  absence  of  the 
transmission  line.  If  the  region  to  the  right  side  of  the  line  in  Eigure  1 1  is  free  space,  then  there 
will  be  no  reflection  of  the  incident  field.  However,  if  there  is  a  local  groundplane  at  some 
distance  from  node  4,  there  can  be  a  reflection  of  the  incident  field,  with  a  corresponding 
reflection  coefficient  being  defined. 

Taking  all  of  the  reflected  terms  into  account,  the  network  reflection  matrix  can  now  be 
written  as 
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or  in  terms  of  the  normalized  voltage  variables, 
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3.3  The  Extended  BET  Equations 

Eqs.(30)  and  (34)  can  now  be  solved  for  the  vector  of  incident  waves  at  the  nodes  using 
exactly  the  same  technique  demonstrated  in  Eqs.(lO)  -  (12).  This  provides  the  following 
extended  BET  equation  for  the  incident  voltages  and  E-fields  at  the  nodes: 
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Using  Eq.(34)  in  the  above  expression  then  provides  the  following  total  voltage  and  E- 
field  extended  BET  equation: 
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4.  Application  of  the  Extended  BLT  Equation  to  a  More  Complex  Problem 

With  the  background  information  on  the  development  of  the  BLT  equation,  and  its 
extension  to  permit  the  inclusion  of  propagating  E-fields  near  the  transmission  line  signal 
paths,  it  is  possible  to  consider  the  more  complex  problem  suggested  in  Figure  10.  Figure  13 
defines  the  various  distances  and  angles  that  enter  into  this  problem.  Note  that  in  this  case,  the 
scattering  bodies  in  the  outside  and  inside  regions  are  includes,  and  the  effects  of  the  scattered 
fields  from  these  bodies  are  incorporated  in  the  BLT  equation. 


Figure  13.  Geometry  of  the  more  general  problem  of  Figure  10,  showing  relevant 
distances  and  angles. 

For  this  more  complex  geometry,  the  signal  flow  diagram  is  shown  in  Figure  14.  In  this 
diagram,  the  various  incident  and  reflected  waves  at  the  nodes  in  the  network  are  illustrated. 
On  the  left  of  the  diagram.  Tube  2  represents  the  propagation  path  of  EM  fields  in  the  region 
outside  the  shield.  These  fields  are  produced  by  an  EM  source,  which  we  represent  here  simply 
as  magnetic  and  electric  dipole  moments.  This  tube  is  connected  to  Node  3,  which  represents 
the  scattering  obstacle  in  the  outside  region,  and  to  Node  4,  which  is  the  aperture  in  the 
conducting  shield.  This  node  acts  much  like  a  node  in  a  transmission  line  network,  in  that  it  is 
represented  by  a  scattering  matrix  that  relates  incident  and  reflected  waves  at  the  junction. 
Unlike  the  transmission  line  junction,  however,  the  waves  in  this  case  are  EM  fields  and  not 
current  or  voltage  waves  on  the  lines. 

As  in  the  simple  model  in  Figure  12,  the  transmission  line  is  represented  by  Tube  1, 
with  end  Nodes  1  and  2.  Tubes  3,  4  and  5  represent  field  propagation  paths,  and  Node  6  is  the 
other  scattering  object  in  the  inside  (shielded)  region.  As  before.  Node  5  is  the  field  coupling 
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node,  in  which  there  is  an  interaction  between  the  incident  fields  on  the  line  and  the  voltage 
traveling  waves. 
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Figure  14.  Signal  flow  diagram  for  the  shielding  problem  of  Figure  10. 

Using  the  same  analysis  approach  as  in  Section  3.3  for  the  field  excited  line,  the 
reflected  voltages  on  Tube  1  at  Nodes  1  and  2  are  given  by  the  following  expression,  which  is 
similar  to  Eq.(25).  Note,  however,  that  in  this  case  there  are  two  field  coupling  terns  involving 
the  incident  E-fields  from  Tubes  3  and  4  and  the  coupling  functions  E^  and  E2.  Observe  that 

these  coupling  terms  are  different  in  that  they  depend  on  the  angles  y/j^  and  y/4,  which  are 
different  as  noted  in  Eigure  13. 
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Eor  Tube  2,  the  relationship  between  the  reflected  and  incident  E-fields  at  nodes  3  and 
4  are  given  by  an  equation  similar  to  Eq.(29),  but  this  time  without  the  transmission  line 
coupling  excitation  terms.  This  equation  is  presented  in  Eq.(44)  and  only  the  excitation  effects 
of  the  primary  sources  are  noted. 
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Tubes  3  and  4  both  provide  the  field  coupling  to  the  transmission  line  (Tube  1)  and 
consequently,  they  will  contain  the  field  coupling  terns  F  again.  The  appropriate  equations  for 
the  reflected  and  scattered  waves  on  these  tubes  are  similar  to  those  of  Eq.(28),  but  without  the 
primary  source  terms,  and  are  given  by 
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for  Tube  3,  and 
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for  Tube  4.  Finally,  the  corresponding  field  propagation  relationships  for  Tube  5  are 
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Equations  42  -  47  can  be  collected  into  a  10x10  propagation  matrix  equation  for  the 
reflected  voltages  and  fields  at  the  nodes,  as  shown  in  Eq.(48).  This  equation  as  the 
counterpart  of  Eq(30b)  for  the  simple  network  of  Figure  12. 
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In  addition  to  the  propagation  relationship  of  Eq.(48),  we  ean  also  develop  a  10x10 
matrix  relationship  arising  from  the  seattering  of  waves  at  the  junctions  -  just  as  done  in 
Eq.(34b)  for  the  simple  problem.  In  developing  this  expression,  however,  we  will  use  a 
slightly  different  notation  from  that  used  in  Eq.(34b).  In  the  following,  we  let  the  superscript 
of  the  reflection  coefficient  denote  the  node  at  which  the  reflection  is  taking  place,  and  the 
subscripts  indicate  the  various  refection  coefficient  at  the  node.  Eor  example,  at  Node  1,  there 
is  only  a  simple  transmission  line  voltage  reflection  coefficient,  and  it  is  denoted  by  pi .  Eor 
this  node,  the  scattering  relationship  is 
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and  the  reflection  coefficient  same  as  that  given  in  Eq.(9).  Similarly,  the  reflection  coefficient 
on  the  transmission  line  at  node  2  is  given  as 
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Node  3  is  also  a  single  branch  node,  but  here  the  E-fields  are  the  important  parameters. 
The  scattering  relationship  is  given  as 


E'i=PlEZ  (51a) 

or  multiplying  by  the  length  parameter  a3  to  put  the  field  variable  into  the  units  of  volts,  we 
have 
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The  reflection  at  Node  4  is  different,  however,  because  there  are  three  tubes  (2,4  and  5) 
connecting  to  the  aperture.  Thus,  we  require  a  3x3  scattering  matrix  of  the  form 
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Details  of  this  scattering  matrix  will  depend  on  the  aperture  size  and  shape,  and  also  depends 
on  the  various  angles  of  incidence  and  reflection  from  the  aperture.  Explicit  expressions  for 
these  terms  remain  to  be  developed  for  simple  apertures,  and  will  be  the  subject  of  a  future 
investigation. 

Node  5  is  the  “field  coupling  node”  and  as  in  the  simple  example  case  it  involves 
reflections  of  the  fields  at  this  node  with  the  transmission  line  removed.  Since  there  are  two 
tubes  (3  and  4)  attached  here,  we  expect  a  2x2  matrix  of  the  form 
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Finally,  node  6  is  a  2-tube  eonneetion  at  the  distant  scattering  body  and  the  pertinent 
expression  for  the  scattering  coefficients  is 
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The  individual  node  scattering  relationships  in  Eqs.(49)  -  (54)  can  be  put  into  a  large 
10x10  matrix  equation.  In  doing  this,  however,  we  must  re-order  the  matrix  entries  so  that  the 
ordering  of  the  incident  and  reflected  voltage/field  vectors  correspond  to  that  in  the 
propagation  matrix  equation  in  Eq.(48).  This  results  in  the  following  system  scattering  matrix: 
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Substituting  this  equation  into  Eq(48)  and  solving  for  the  incident  voltage  vector  gives 
the  first  form  of  the  BLT  equation  for  the  incident  response  components,  as  reported  in 
Eq.(56).  Adding  Eqs.(55)  and  (56)  yields  the  total  responses  at  the  nodes,  as  given  by  Eq.(57). 
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5.  Summary 

This  paper  has  outlined  an  analysis  methodology  for  treating  EM  scattering  and 
aperture  penetration  problems  with  the  formalism  of  the  BET  equation.  In  addition  to  the  usual 
transmission  line  paths  on  which  voltage  and  current  waves  propagate  with  little  attenuation,  in 
the  extended  problem  has  EM  “signal  paths”  that  have  a  1/r  attenuation,  and  which  are  able  to 
couple  to  transmission  line  structures  along  the  length  of  the  line,  not  just  at  the  end  nodes. 

To  analyze  this  field  coupling  phenomenon,  we  introduce  the  concept  of  a  “field 
coupling  node”,  which  takes  into  account  the  incident  EM  field  coupling  to  the  transmission 
line.  Moreover,  by  reciprocity,  this  node  also  provides  effects  of  EM  field  radiation  by  the 
transmission  line  back  to  the  source  or  other  scatterers  in  the  problem  space. 

To  illustrate  this  methodology,  two  examples  have  been  given:  a  simple  radiating 
source  coupling  to  an  isolated  transmission  line,  and  a  source  near  an  aperture  in  a  shield, 
coupling  to  a  line,  with  several  scattering  objects  nearby.  In  both  of  these  cases,  only  the 
analysis  methodology  has  been  outlines;  numerical  calculations  and  verification  of  this 
analysis  have  yet  to  be  done,  and  will  be  the  subject  of  a  subsequent  note. 
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Appendix  A 

Field  Excitation  of  a  Transmission  Line  using  the  BLT  Equation 


In  this  appendix  we  wish  to  derive  the  BLT  equations  for  the  load  voltages  and  currents 
for  the  case  of  distributed  field  excitation  of  the  transmission  line,  as  portrayed  in  Figure  Al. 
As  discussed  in  [Al],  this  excitation  can  be  represented  by  the  tangential  incident  E-field  along 
the  line  and  at  the  end  conductors  of  the  terminating  impedances.  Generally,  the  distributed 
voltage  sources  will  induce  both  antenna  mode  and  transmission  line  mode  responses,  but  as 
previously  mentioned,  we  will  consider  only  the  transmission  line  responses  in  this  discussion. 


Figure  Al.Field  excitation  of  a  two-wire  transmission  line,  showing  the  voltage  excitation 
sources  needed  in  the  BLT  equation. 

The  transmission  line  excitation  along  the  length  of  the  line  can  be  represented  by  a 
distributed  voltage  source,  V'(x)  which  is  given  by  the  difference  of  the  v-components  of  the 
incident  E-field  along  the  line.  This  is  expressed  as 


VXx)  =  E‘:\x,d)-E‘:^(x,0).  (Al) 

In  addition,  there  are  two  lumped  voltage  sources  at  either  ends  of  the  line,  and  Vg2, 
which  take  into  account  the  interaction  of  the  incident  E-field  with  the  vertical  conductors  in 
the  loads  at  the  ends  of  the  line.  These  are  given  in  terms  of  the  z-components  of  the  incident 
E-field,  evaluated  at  z  =  0,  and  are  expressed  as 

-  Ef  (0, 0)d  and  -  Ef  (E,  0)d  .  (A2) 

The  load  responses  given  by  the  BET  equations  (14)  or  (16)  in  the  main  body  of  this 
report  can  be  thought  of  as  Green’s  functions  for  computing  the  responses  due  to  a  distributed 
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field  excitation.  Setting  the  current  excitation  term  to  zero  in  Eq.(14),  the  following  Green’s 
function  for  the  load  voltages  due  to  a  lumped  voltage  source  at  location  v  can  be  defined: 


G,  (xh 

1  + A 

0 

.  “A  1 

-1 

Gj  (x)\ 

0 

l  +  Pi 

-P2I 

(A3) 


Here  Gj{x)  denotes  the  Green’s  function  for  the  load  voltage  response  at  node  1  and  G2(x)  is 
the  Green’s  function  for  the  voltage  at  node  2.  Similar  Green’s  functions  for  the  load  currents 
can  be  defined  using  Eq.(16). 

It  is  clear  that  by  inserting  in  the  voltage  source  distributions  given  by  Eqs.(Al)  and 
(A2)  into  Eq.(A3)  and  integrating  over  the  transmission  line,  the  only  modification  to  the  BET 
equation  is  a  change  of  the  source  terms.  In  this  manner,  the  BET  voltage  equation  for 
distributed  source  excitation  becomes 


>1' 

'l  +  A 

0 

-Pi 

-1 

_^2: 

0 

1+ A_ 

1 - 

1 

to 

1 _ _ 

(A4) 


where  the  source  vector  is  now  given  as 


-  Ef  (x,0))e’*'  dx  +  Ef  (0,0)  J  -  Ef  (E,0)  J 
Ef  (v,0)y^^“"‘W-Ef  (0,0)J  +  Ef  (E,0)J 


.  (A5) 


Note  that  in  this  equation,  the  first  term  in  the  source  vector  is  the  integral  term  accounting  for 
the  effects  of  the  distributed  field  excitation  of  the  line,  and  the  second  and  third  lumped  terms 
account  for  the  excitation  of  the  end  conductors  of  the  line. 

Eor  the  special  case  of  a  vertically  polarized  plane  wave  of  amplitude  E'”*^  propagating 
in  the  x  -  z  plane  with  an  incidence  angle  if/  (as  shown  in  Eigure  Al),  the  x  and  z  components 
of  the  incident  E-field  can  be  written  as 


Ef  (x,  z)  =  sin  _ 

(Aba) 

Ef  (v,  z)  =  E'”"  cos 

(Abb) 

Assuming  that  the  transmission  line  propagation  constant  is  the  same  as  that  of  free 
space  (e.g.,  y=  jk),  inserting  the  above  expressions  for  the  incident  E-fields  into  Eq.(A5), 
performing  the  indicated  integrations  and  then  combining  terms  provides  the  following 
expression  for  the  source  vector  for  plane  wave  excitation  of  the  line: 
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^^jkL(l-cosyf) 

2 

^ikL^-jkUl+cosyr)  _ 

=  E 


inc 


F,m, 


(A7) 


This  last  equation  serves  to  define  the  field  coupling  functions  Fj  and  F2,  which  define  how 
strongly  the  traveling  waves  are  excited  by  the  incident  field.  Note  that  ref.[l]  presents  the 
more  general  expression  for  this  source  vector  for  the  case  of  an  arbitrarily  incident  and 
polarized  EM  field  acting  on  the  line  and  for  the  case  of  the  transmission  line  propagation 
constant  ^'differing  from  the  free  space  propagation  constant. 

As  in  the  case  of  a  single  lumped  excitation  source  on  the  line,  the  Sj  source  term  in 
Eq.(A7)  can  be  attributed  to  the  positive  traveling  wave  on  the  line,  as  excited  by  the  incident 
field,  and  the  S2  term  is  due  to  the  excitation  of  the  negative  traveling  wave. 
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Appendix  B 


Radiation  Fields  Produced  by  Traveling  Wave  Components  on  a 

Transmission  Line 


One  of  the  elements  needed  in  the  generalization  of  the  BLT  equation  to  include  field 
propagation  effects  is  the  radiation  produced  by  traveling  waves  on  a  transmission  line.  As 
noted  in  [Bl],  this  can  be  determined  using  the  concepts  of  electromagnetic  reciprocity. 
Alternatively,  it  is  possible  to  simply  integrate  over  the  current  distribution  on  the  line,  which 
arises  from  the  traveling  wave  currents.  This  latter  approach  is  considered  in  this  appendix. 

Figure  Bl  shows  a  transmission  line  of  length  L  with  a  positive  traveling  voltage  wave 
y+(v).  Not  shown  in  this  diagram  is  a  corresponding  negative  traveling  wave.  Both  of  these 
waves  are  excited  in  some  manner  by  some  sort  of  an  excitation  function,  and  the  behavior  of 
these  waves  is  determined  by  the  termination  impedances  and  the  excitation  function.  At  this 
stage,  we  are  not  interested  in  the  complete  solution  of  the  voltage  or  current  on  the  line  (as 
this  ultimately  arises  from  a  solution  of  the  BLT  equation),  but  we  are  seeking  to  determine  the 
radiated  field  contributions  from  each  of  these  traveling  waves  -  expressed  in  terms  of  the 
incident  voltage  waves  at  each  end  of  the  line. 

Radiated  field 


>2  Z\ 


ox  L 

Node  #1  Node  #2 


Figure  Bl.  Illustration  of  a  positive  traveling  wave  on  a  transmission  line,  with  the 
production  of  a  radiated  field  . 

The  radiated  E-field  is  assumed  to  be  in  the  direction  at  a  distance  from  the  origin 
of  the  coordinate  system,  as  shown  in  the  figure.  We  assume  that  is  sufficiently  far  from  the 
line  so  that  the  radiated  field  is  a  plane  wave,  with  primary  component  E^. 

For  the  positive  traveling  voltage  wave  y+(v),  the  corresponding  current  traveling  wave 
is  given  as  I+(x)  =  V+(x)IZ^,  where  Z^.  is  the  characteristic  impedance  of  the  transmission  line. 
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For  this  traveling  current  wave,  there  are  individual  traveling  wave  currents  on  each  conductor 
of  the  line  given  by  Ij  and  I 2,  as  shown  in  Figure  Bl,  with  the  constraint  that  I 2  =  -I]- 

To  determine  the  E- field  radiated  by  the  traveling  wave  currents  on  this  line,  it  is 
necessary  to  recall  the  fundamental  radiation  formula  for  the  field  produced  by  an  infinitesimal 
current  element  Idl,  as  shown  in  Figure  B2.  In  the  far-zone,  the  radiated  E-field  is  given  as 
[Bl] 


Eo  = 


4;rr„ 


sin  0  Idl  e 


(Bl) 


z 


I  dl 


r 


0 


Figure  B2.  Radiation  field  geometry  for  a  point  current  element,  /  dl. 


Noting  that  E^is  in  the  opposite  direction  than  the  field  component  in  Eigure  Bl, 
and  that  sin  0  =  sintj/,  the  differential  radiated  E-field  component  due  to  the  pair  of  currents  at 
location  x  on  the  line  may  be  expressed  as 


dE,  =  -  sin  w 
^  Ak 


p-M  -jkr^ 

I^(x) - +  l2(x) - 

h  ^2  y 


(B2) 


Using  the  fact  that  the  traveling  wave  currents  can  be  expressed  in  terms  of  the  incident 
traveling  voltage  wave  at  node  2  (at  x  =  L)  as 


/i(v)  =  -/2(v) 


^-r{x-L) 


together  with  the  far-field  approximations 


1  1  1 


jk{r^+xco^\(f-d  ^\x\y) 


—  ^-jk{r„+xcosl//) 


(B3) 


(B4a) 

(B4b) 

(B4c) 
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the  differential  field  becomes 


•  1  T  7  inc 

dE,  =  -l).  (B5) 

^  Am-,  ^ 

Noting  that  -1)=  jkdsinl//  and  assuming  again  that  ;^=  7^,  Eq.(B5)  can  be  put 

into  the  following  form: 

h-^  r!  7  ~ 

dE,  =—^- - sin’  at  .  (B6) 

^  4jrZ^  r„  ’ 


To  obtain  the  total  contribution  to  the  radiated  E-field  from  the  traveling  wave  currents, 
Eq.(B6)  must  be  integrated  over  the  length  of  the  line.  This  is  expressed  as 


E\=^  E^{x)dx  = 
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(B7) 


The  radiated  field  in  Eq.(B7)  is  not  the  complete  field  from  the  positive  traveling  wave, 
however,  because  the  contributions  of  the  traveling  wave  components  in  the  end  loads  have  not 
been  included.  As  shown  in  Eigure  B3  there  is  a  current  Ij  flowing  in  the  end  of  the  line  at  v  = 

0,  and  a  current  I2  at  the  end  at  x  =  L.  These  currents  may  both  be  expressed  in  terms  of  the 
incident  voltage  wave  at  node  2  as 


and 


(B8a) 
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(B8b) 
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Figure  B3.  The  traveling  wave  currents  at  the  vertical  ends  of  the  transmission  line. 

Using  the  expression  (Bl)  again,  the  E- field  contribution  from  these  end  currents  can 
be  written  as 


- 
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](mo 

imo 


V  /  \ 

~j^fn  '2  ikLi-i  - jk{l+cosw)L  \ 

cosi//e  - e  \l  —  e  ^  ). 
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(B9) 


Combining  Eqs.(B7)  and  (B9)  provides  the  total  radiated  field  from  the  positive 
traveling  wave  component  as 


-jkr„ 
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Using  the  definition  of  the  backward  wave  field  coupling  function  F2,  as  defined  by  Eq.(A7)  in 
Appendix  A,  the  radiated  field  for  the  positive  traveling  wave  on  the  line  given  in  Eq.(BlO)  can 
be  written  as 


'h  7  ~J^^o 

El  - F.iipwr  ,  (BIO) 

^  In  Z^  ^ 

which  is  essentially  a  statement  of  reciprocity  between  the  reception  and  transmission 
properties  of  the  transmission  line. 

In  a  similar  manner  to  the  development  above,  the  E-field  produced  by  the  traveling 
wave  components  on  the  line  and  in  the  end  loads  can  be  computed.  In  this  instance,  the 
expression  will  be  developed  in  terms  of  the  incident  voltage  wave  at  node  1,  yj'"'"to  ultimately 
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be  compatible  with  the  requirements  of  the  BLT  equation.  Upon  carrying  out  the  required 
integration  and  simplification  of  results,  the  E-field  from  this  negative  traveling  wave 
(including  the  end  current  effects)  is 
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Vine  (  ji 

1 


Jk{l  -cosy/)L 


-i) 


sin^  Xj/ 

(l  -  cos  y/) 


-cosy/ 


Ak  Z, 
jk  Z,  c' 
l/t  Z„ 


ro 

■Jkr„ 


y  inc  ik(\-cos  y)L  _  j 


(Bll) 


vrF,{y/) 


where  F](y/)  is  the  field  coupling  function  for  positive  traveling  waves  on  the  receiving 
transmission  line,  given  by  Eq.(A7)  of  Appendix  A. 

By  combining  Eqs.(BlO)  and  (Bll),  the  total  E-field  radiated  by  the  two  traveling 
waves  on  the  line  can  be  determined: 

E^=E^  +  E-=^^—  {vrF,(iE)  +  VrF,(y/)).  (Bll) 

In 


References 


BE  Tesche,  E.  M.,  et.  ah,  EMC  Analysis  Methods  and  Computational  Models,  John 
Wiley  and  Sons,  New  York,  1997. 


43 


